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We estimate the intrinsic charm contributions to the quark flavor and spin observables of the 
nucleon in the SU(4) quark meson fluctuation model. In this model, the charm or anticharm 
reside in the charmed mesons created by the nonperturbative quantum quark-meson fluctuations. 
The intrinsic charm content in the proton, 2c/ Y^{q + q) ~ 0.011 ± 0.008, is almost one order of 
magnitude smaller than the intrinsic strange content. The intrinsic charm helicity is also small and 
negative, Ac ~ —(0.009 ± 0.006). The fraction of the total quark helicity carried by the charm is 
|Ac/AEj ~ 0.021 ±0.014. The ratio of the charm with positive helicity to that with negative helicity 
is C|/c| = 35/67. For the intrinsic strange component, one has Sj/sj ~ 7/13. A detail comparison 
of our predictions with data and other models or analyses is given. The intrinsic charm contribution 
to the EUis-Jaffe sum rule is also discussed. 



PACS numbers: 14.65.Dw,14.20.Dh,12.39.Fe,11.30.Hv 

1. Introduction 

As suggested by many authors long time ago |j, 
there are so called 'intrinsic ' heavy quark components in 
the proton wave function. The intrinsic heavy quarks 
are created from the quantum fluctuations associated 
with the bound state hadron dynamics. They exist in 
the hadron over a long time independent of any exter- 
nal probe momentum transfered in the collision. The 
probability of finding the intrinsic heavy quarks in the 
hadron is completely determined by nonperturbative 
mechanisms. On the other hand, the extrinsic heavy 
quarks are created on a short time scale in association 
with a large transverse momentum reaction and their 
distributions can be derived from QCD bremsstrahlung 
and pair production processes, which lead to standard 
QCD evolution. The study of nonperturbative mecha- 
nism <75 — > 99 shows that the intrinsic heavy quark 
contribution scales as l/mg, where mg is the mass of 
the heavy quark, the probability of finding the intrinsic 
bottom is expected to be much smaller, hence we will 
only consider the intrinsic charm (IC) component, e.g. 
\uudcc >, in the nucleon. 

Many works on intrinsic charm are based on models 
of the nucleon (see e.g. MIT bag model NJL model 
1^ , meson cloud model [|l|, ^ , instanton model ^ and 
others 0, etc.) or combination of the IC hypothesis and 
analysing the existing DIS data on charm production to 
obtain some information of the intrinsic charm content 
(see e.g. [| | Jp, |ll|, |l|, [l| ^ 0, |l|). As 
pointed out in fl^ that a non-perturbative IC component 
was unable to explain enhancement in the e+p neutral 
current cross section at HERA kinematics, but HERA 



data may provide a good test of the intrinsic charm if 
the charm component of F2 is directly isolated and no 
new physics exists. Another discussion given in p6) sug- 
gested a nonsymmetric charm distribution, for which c is 
much harder than c, which may account for the HERA 
anomaly, a significant enhancement of cross sections at 
HERA at large x and However, the most 

recent analysis given in U7L which includes both the 
perturbative photon-gluon fusion process and nonpertur- 
bative intrinsic charm contribution, found no conclusive 
evidence for the intrinsic charm and only provided an up- 
per bound of around 0.4%. Therefore what is the effect 
of charm in the structure of the proton is still an open 
and interesting question. A comprehensive review on the 
nucleon sea, which includes light quark, intrinsic strange 
and charm, can be found in [ pT| . 

The chiral quark model or more precisely the quark 
meson fluctuation model (in SU(2) version) was first 
suggested in to explain the sea flavor asymmetry, 
d — u > 0. The model was extended to SU(3) with- 
out and with SU(3) symmetry breaking to describe the 
quark spin structure (see e.g. p^, pil 
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references therein), and orbital structure (see e.g. |2|] 
and references therein) of the nucleon. In this model, the 
nucleon structure is determined by its zeroth order va- 
lence quark configuration and all possible quantum fluc- 
tuations of the valence quark into a recoil quark and a 
meson. It provides a natural and nonperturbative expla- 
nation of a nonzero (intrinsic) strange sea and negative 
strange sea polarization through the quark-kaon fluctua- 
tion (e.g. u — > s + K'^). The results agree quite well with 
the data obtained in deep-inelastic scattering and other 
related experiments. 

Based on earlier encouraging result we suggested 
an SU(4) version of quark-meson fluctuation model and 
predicted a nonzero intrinsic charm in the proton [ p9[ . 
In the SU(4) version, charm or anti-charm quarks re- 
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side in the charmed mesons which are created by non- 
perturbative quantum quark-meson fluctuations such as 
u ^ c + . These charm or anti-charm quarks are es- 
sentially intrinsic components of the hadron. 

The main motivations of extending to SU(4) version 
are: 1) There is no physical reason to forbid the fluctua- 
tion from a light quark to a charmed quark and charmed 
meson. Hence it is interesting to investigate what are 
the changes of the SU(3) predictions after including the 
charm or anticharm contributions. 2) According to the 
symmetric GIM model |Q, one should deal with the 
weak axial current in the framework of SU(4) symme- 
try. It implies that the charm quark should play some 
role in determining the spin, flavor and orbital structure 
of the nucleon. Therefore it is natural and reasonable 
to extend the SU(3) version to SU(4) by including the 
charm and anticharm contributions. 

In this paper, we will discuss the effect of the IC con- 
tribution to the flavor and spin observables of the pro- 
ton in the SU(4) quark meson fluctuation model with 
symmetry-breaking. The analytic and numerical results 
of the flavor and spin contents for each quark flavor u, d, 
s, and c in the proton are presented. Predictions are com- 
pared with the existing data and other model results or 
analyses. The intrinsic charm contribution to the EUis- 
Jaffe sum rule is also discussed. 
2. SU(4) model with symmetry breaking 

In the framework of SU(4) quark model, there are six- 
teen pseudoscalar mesons, a 15-plet and a singlet. In this 
paper, the contribution of the SU(4) singlet will be ne- 
glected. It is easy to show that the effective Lagrangian 
describing interaction between quarks and the mesons in 
the SU(4) case is 



( 



where Z?+ — {cd), D 
= (cs), and Dj = 
nents G°i^^^ and G^^ 



(cd), 




= (cu), 

= (cs). The neutral charge compo- 
are defined as 



(2) 

G°^-if^/2^3 + v"VCV^-V°V^ (3) 



(4) 



with 



n° ^ {uu-dd)/V2, ^ {uu + dd - 2ss)V6, (5) 



T]'" = (mm + dd + ss) /VS, ril = (cc) . 



(6) 



Similar to the SU(3) case, we define a = jgisp, which de- 
notes the transition probability of splitting u{d) — s- d{u) + 



7r+( \ Hence ea, e^a, C'^a, and CcO denote the probabil- 
ities of splittings u{d) ^ s + K~^^\ u{d) u{d) + rf , 
u{d) u{d)+i]"^ and u{d) c+D^{D^) respectively. If 
the splitting probability is dominated by the mass effects, 
we expect < < C'^a < ea ~ e,,a < a, or 



< ec < C < e 



< 1. 



(7) 



It implies that the probability of emitting heavier me- 
son such as D from light quarks is smaller than that of 
emitting the lighter mesons such as K, t], and 77', etc. 
Since the charm contributions only appeared in the ob- 
servables which are explicitly related to the charm or 
anticharm variables, therefore we expect that the SU(4) 
model will not significantly change the SU(3) predictions 
on observables which are not explicitly charm-related. 
On the other hand, the SU(4) version leads to many new 
predictions on observables which are directly related to 
the charm or anticharm. These predictions can be tested 
by the future experiments. Before going to further dis- 
cussion, we would like to make a few remarks. 

• The quark contents of the ry" and r/"^ shown in 
Eqs.(5) and (6) imply that they are actually rjs and 
771. Hence the octet-singlet mixing is neglected in 
our description. 

• We have shown in ||2^ that a better SU(3) descrip- 
tion can be achieved by taking e,, ~ e and C'^ << 1, 
i.e. the effective U(l)-breaking parameter is smaller 
than other SU(3)-breaking parameters. This result 
was deduced from the data of d — m, which requires 
C' should be small and negative. More detail dis- 
cussion can be found in section II and Fig.l in [p7| . 
It can be shown that including the charm contri- 
butions will not change this conclusion. Therefore, 
we may choose C'^ = as one possible option in 
the numerical calculation, i.e. neglect the singlet 
contribution from the beginning. In this case, the 
number of independent parameters would be only 
three: a, e, and Ec- Since our analytical formulae 
are not restricted to this approximation, we can 
also discuss C'^ 7^ case (see discussion in section 
4 below). 

• It should be noted that the definition of G", in 

u{a) 

SU(4) case is different from that in the SU(3) case. 
Besides an additional charm term, —rj'^^/e^JW, the 
coefficient of 77° term is changed. For G°, even the 
coefficient of 77'° term is also changed. Therefore, 
the SU(4) formahsm cannot be reduced to SU(3) 
simply by taking e,, only [see later discussion 
on Eqs. (15) and (16)]. The G° is completely new 
in the SU(4) version. 

• As pointed out in the original chiral quark model 
[H H m m, H, the nucleon properties 
are defined in the scale range between Aqcd {'^ 
0.2 - 0.3 GeV) and A^sb (~ 1 GeV), where the 
spontaneous breaking of chiral symmetry leads to 
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the existence of Goldstone bosons. In the quark 
meson fluctuation model (SU(3) in [|8) and SU(4) 
in this paper), we will assume all calculated quan- 
tities are also defined in the same scale range [(0.2 
GeV)2 < fj.^ <(1.0 GcV)2]. At this range, the 
sea content should be dominated by the intrinsic 
component created by nonperturbative mechanism, 
quark-meson fluctuation, |3l[| and the extrinsic sea 
component is expected to be small. 

• We note that in our formalism, only the integrated 
flavor content [q{Q^) — dxq{x,Q^)] and helicity 

content [Aq{Q'^) — dxAqlx^Q"^)] are discussed. 
To make brief of the formalism, we will omit the 
argument everywhere unless specified otherwise. 

In addition to the allowed fluctuations discussed in the 
SU(3) valence quark is now allowed to split up or 

fluctuate to a recoil charm quark and a charmed meson. 
For example, a valence u-quark with spin-up, the allowed 
fluctuations are 

^ (9) 

W| (10) 

Similarly, one can list the allowed fluctuations for 
d-f, dj, Sf, and S|. Since we are only interested in 
the spin-flavor structure of the nucleon (or other non- 
charmed hadrons), which does not have valence charm 
quark, hence the fluctuations from a valence charm quark 
will not be discussed. 

The spin-up and spin-down quark or antiquark con- 
tents in the proton, up to first order of the quantum 
fluctuation, can be written as (a detail SU(3) version see 
e.g. H), 

np(g^_^,or g^_^) = "p°^(«'*)-^9"('?Ta'°'^ 

q=u.d h=1 ,1 

(11) 

where q' = u, d, s, c and n^\qi^i) are determined by the 
zeroth order, i.e. naive quark model (NQM), valence 
quark wave function of the proton and 

(ut) = 5/3, n(")(«J = l/3, 

(12) 

<)(dT) = 1/3, n(")(di) = 2/3. 

In Eq.(ll), Pq^ ^{q'-^ ^) and Pq^^iiq^^i) are probabilities 
of finding a quark q'^ ^ or an antiquark q'^ ^ arise from 
all quantum fluctuations of a valence quark qi^i- The 
probabilities Pq^ ^{q'^ ^) and Pq^ ^(^t obtained 
from the effective Lagrangian (1) and have been listed 
in Table I, where only Pq^ (q'^ ^) and Pq^ (q'^ ^) are shown. 
Those arise from q^ can be obtained by using Pq^ {q\ = 



Pq^ (^j. I ) and Pq^ (q'] i) = Pq^ (^J | ) . The notations ap- 
peared in Table I are defined as 

/ = l/2 + e^/6 + C'V48 + e,/16, 

(13) 

= 2e^/3 + C'V48 + ec/16, 

and 

i=l/2 - ^/6-C7l2, S = -^/3 + C7l2, 

(14) 

C EE 2^1-5 + C' 11^, D=^,/A. 

Analogous to the SU(3) case, the special combinations 
A, B, C, and D stem from the definitions (2)-(4), which 
show the quark and antiquark contents in the neutral 
bosons G^^^ ,.^^. The numbers fa and /^a stand for the 
probabilities of the quark splitting M-f(£i-f) "1(^1) + 
G°^^^ and S| -I- respectively. 

To reduce to the SU(3) version, one should take the 
limit ec ^ and change C' to 4^' simultaneously. This 
is because we have different neutral combinations G° ^ ^ 
from the SU(3) ones. This difference shows that the 
charm contribution is not only presented in the charm 
term which is proportional to ec, but also affected by 
the C term. Taking Cc ^ and C' 4^', one can see 
from Eq. (13) that the / and fs indeed reduce to the 
corresponding quantities in the SU(3) case, 

if)sui3) = l/2 + e„/6 + C'V3, 

(15) 

ifs)sui3) = 2e„/3-fC'V3, 

which is exactly Eq.(6b) in Taking ec ^ and 

C 4C', one has I) —> and 

A^Asu(3)/S, B^Bsui3)/i, C^Csui3)/i, (16) 

where Asu(^^), Bsm^^, and £"5(7(3) are the same as A, B, 
and C in Eq. (6a) in [||. 

If Eq.(12) is replaced by ni^Vt) = 5/3, nl°\di) = 
1/3, nn\u^) = 1/3, and nn\ui) — 2/3, from (11), we 
will obtain all flavor and spin contents in the neutron. 
3. Flavor and spin contents in the proton 

We note that the quark flips its spin in the quark split- 
ting processes (ZT,(i) ~^ processes in (8) 
and (9), but not in the (10) (no splitting). The quark 
helicity non-flip contribution in splitting processes is en- 
tirely neglected. This is a basic assumption in the model 
and seems to be consistent with the picture given by the 
instanton model. 

3.1. Quark flavor content 

Using (11), (12) and the probabilities Pq^ ^{q'^ ^) and 
Pq^ i i) listed in Table I, one obtains all quark and 
antiquark flavor contents in the proton, 

u^2 + u, d^l + d, s = + s, c = + c, (17) 
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where 

u = a[l + A^+ 2(1 - Af], a[2(l + A^) + (1 - A)\ 

(18) 

s = Za[€ + B% c = 3a[ec + i)2]. (19) 
From (18), one obtains 

M/(i= 1 - 6i/[(3i- 1)2 + 8], (20) 

and 

d-u = 2aA. (21) 

Similarly, one can obtain 2c/ (u+d), 2c/ (u+d), 2c/ 
q) and other flavor observables related to charm quark 
and cahrm antiquark. We note that in the SU(4) case, 
one has 

^(g + g-) = 3(l + 2aa), (22) 

9,9 

where = 1 + e + + /• It is easy to verify that in the 
limit Ec 0, and using Eq.(16), all SU(4) results reduce 
to those given in the SU(3) case. Two remarks on (20) 
and (21) are in order. 

• Defining the ratio, r = u/d, we obtain from Eq. 
(20) that A = 2-r±Vl4r - Sr"^ ~ 5. Since A must 
be a real number due to (21), hence 14r — 'Hr^ — 5 
must be positive or zero, and the ratio r satisfies 

1/2 < u/d < 5/4, (23) 

which seems to be consistent with the experimental 
data shown in Tabic 111. 

• Since |C'/2| < ^/Cjj in our approximation, from (14), 
we have ^ < 1/2 (note that C is negative) and 

d-u<a. (24) 

It gives a lower bound of the parameter a — the 
probability of quark splitting q q' + tt. 

3.2. Quark helicity content 

Similar to the SU(3) version, we have 

q'—u,d 

(25) 

Using Eqs. (12), (25), and Table I, we obtain 

Au = 4[l-a(e + ec + 2/)]/3-a, (26) 

Ad = -[l-a(e + ec + 2/)]/3-a, (27) 

As = -ae, (28) 



Ac = -aec, (29) 
AI]= Ag = l-2a(l + e + £, + /), (30) 

q—u,d^s,c 

or 

AS/2 = 1/2 (31) 

where is defined in (22), and 

Ag = 0, (g u,d,s,c). (32) 

Before going to the numerical calculation, we would like 
to make several remarks on some results which are basi- 
cally independent of parameters. 

• Comparing with the SU(3) case, a new Cc term, 
which presents the intrinsic charm contribution, ap- 
peared in Am, Ad, Ac, and AS, but not in As. 
This is because there is no process which can mix 
the strange helicity and charm helicity contribu- 
tions. As can only come from the processes like 
U| — !■ S| -f , while Ac comes only from the pro- 
cesses like ^ ci + . Although there are s, s 
(or c, c) in all neutral bosons G° d s they give no 
contributions to As (or Ac) due to S| and S| (or 
C| and C|) appeared with equal probability in these 
neutral bosons. 

• The charm quark helicity Ac, (29), is nonzero as far 
as Ec is nonzero. Analogous to the strange quark 
helicity, Ac is definitely negative^ but the size of the 
intrinsic charm helicity depends on the parameter 
ec and the splitting probability a. 

• The physical meaning of (31) is that the total loss 
of the quark helicity arises from four splitting pro- 
cesses with quark spin-flip, three in (8) and a new 
splitting in (9). Comparing with the SU(3) case, 
where AE/2 = 1/2 — a(l -I- e -I- /), we now have an 
additional reduction, — aec, of the total quark spin 
due to the splitting related to the charm. 

• In the splitting process u^{i) C|(|) -I- Z)", the 
anticharm resides only in the charmed meson, e.g. 
D^{c, u). The probabilities of finding Cf and C| are 
equal in the spinless charmed meson. Therefore 
Ac — 0. Similar discussion in the SU(3) case pi] 
has led to Aq — for q — u, d, s. The result (32) 
shows that the helicities of the sea quark and anti- 
quark are not equal, Aqsea 7^ Aq. This is different 
from the usual gluon splitting {g ^ q + q) model 
and gg ^ q + q model |jl|. In the gluon splitting 
and gluon fusion models, the sea quark and anti- 
quark with the same flavor are perturbatively or 
nonperturbatively created as a pair from the gluon 
or gluons and Aqsea = Aq. The DIS data |4^ 
seems to support the prediction Aq ~ but with 
large errors. 
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From (19) and (29), using D"^ = 6^/16 in (14), one 
can see that the ratio 



Ac/c= -16/51 



(33) 



is a constant independent of any sphtting param- 
eters. This is a special prediction for the charm 
flavor in the SU(4) quark meson model. 

• For the strange flavor, from B in (14) and s in (19), 
one has s = s ~ ea(10/3)[l + |C'|/(20VF)] due to 
C'^ << 1. In the Umit C'^ ^ 0, we have 

As/s^-3/10. (34) 

Hence As/s is also a constant in this hmit. 

• From (33) and (32), one obtains the ratio of the 
charm with positive hehcity to that with negative 
hehcity is also a constant. 



C|/q = 35/67- 0.522. 



(35) 



It shows that more is created than C| in the 



splitting processes 



D° and 



-i(T)" 



D . This is because the total con- 



tribution from M| and d| splittings is larger than 
that from and d^ splittings. Similar situation 
occurs for the strange quarks. From (34) and (32), 
we have, for (' = 0, 



s^/si = 7/13 - 0.538. 



(36) 



• For the u-flavor and d-flavor, we do not have similar 
exact relations like (33) and (34). This is because 
the quantities U| — U| and d^ — di depend on t^, 
while U| +wx and d^ + rfx not. Hence Au/u and 
Ad/d depend on all three splitting parameters. The 
numerical calculation (see next section) shows that 

Au/u~ 0.383, Ad/d ~ -0.287. (37) 

From (37) and (32), we obtain 

d|/dx~ 0.554, (38) 

and 

u^/ui ~ 2.241. (39) 

Comparing with the zeroth approximation (NQM), 
where d^/di = 1/2 and u^/ui — 5/1, hence the 
quark-meson splittings lead to a small enhancement 
for d^/di, but a signiflcant reduction for u^/ui. 

• The ratios of the flavor or spin contents predicted in 
this paper are also defined at = fi^ scale, where 
the sea (light or heavy) quark contents are 'in- 
trinsic' and determined by nonperturbative mech- 
anism. To compare the predictions with data at 
higher range, one needs to use perturbative 



QCD and treat these ratios (at = fi^) as in- 
puts, i.e. boundary conditions, in the QCD evolu- 
tion equations. If the numerator and denominator 
of the ratio have the same or almost the same factor 
of dependence, then this ratio will not sensitive 
to the change of Q^. 

3.3. Ellis-Jaffe sum rule 

In the framework of SU(4) parton model, the first mo- 
ment of the spin structure function g^(x, Q^) in the pro- 
ton is 

V\ ^ [\p{x,Q')dx 
Jo 

= (4Au + Ad + As-h4Ac)/18, (40) 

which can be rewritten as 
»i 







giix,Q )dx 



(41) 



= (3a3 -I- as - ai5 + 5ao)/36 
where the notations 

03 = Au - Ad, as = Au + Ad - 2As (42) 



ai5 



Au + Ad + As - 3 Ac 



ao^ Au + Ad + As + Ac 



(43) 



(44) 



have been introduced. Using exchange Au Ad, or 
as ^ — fla in (41), we can obtain F" for the neutron. 
4. Numerical results and discussion. 

Since the probability of charm-related splitting (9) is 
much smaller than u-,d-,s-related splittings (8), we ex- 
pect that the values of parameters a and e in SU(4) 
should be very close to those used in the SU(3) version. 
We choose a = 0.143, e = 0.454 [note that a = 0.145, e = 
0.460 in SU(3)], and leave Cc as a variable, then express 
the quark fiavor and helicity contents as functions of Ec. 
To determine the range of , we plot a^ as the function of 
Ec in Fig.l. Using the most precise data from the neutron 



/?-decay |33 
obtain e, 



ag = [UA/i^v)n^p = 1.2670 ± 0.0035, we 
0.06±0.02. To include other possible theoret- 
ical uncertainties arise from the model approximations, 
however, we prefer to introduce a larger uncertainty and 
take 



0.06 ±0.04. 



(45) 



In addition to the three-parameter set {a, e, ed 
[SU^''-'(4)], we also consider the four-parameter set: 
{a,e,C',ec} [SU('')(4)] to show the 77'" effect. The pa- 
rameter sets used in this paper and in previous SU(3) 
version are listed in Table II. 

Using the parameter sets SU('')(4) and SU(^)(4), the 
flavor and spin observables in the proton are calculated 
and listed in Tables HI and IV respectively. For compar- 
ison, we also list the existing data, the results from pre- 
vious SU(3) description p7|, p8[, the naive quark model 
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(NQM) prediction, and results given by other models or 
analyses. Including the QCD radiative corrections, the 
first moments of the spin structure functions and 5" 
are also shown in Table IV. We also calculated two re- 
duced matrix elements F and D, and the ratios Ga/Gv 
of the hyperon beta decays. The results are listed in Ta- 
ble V. One can see that both SU^'^) (4) and SU^''^ (4) ver- 
sions satisfactorily describe almost all the existing data. 
Furthermore, they give many new predictions on charm- 
related observables, which are in the bold type shown 
in Tables III and IV (for simplicity, we use Cc = 0.06 
in SU^''^(4) without the uncertainty). A few main pre- 
dictions selected from SU''*^(4) results are listed in the 
following: 

2c/ ^(g g) ~ 0.011 ±0.008 (46) 
measures the size of the intrinsic charm in the proton, 



and — t(z) 



ilz/{\ - z). In Eq.(50), 



Ac - -0.009 ±0.006 
is the charm helicity, and 

Ac/as ~ -(0.021 ±0.014) 



(47) 



(48) 



is the fraction of total quark helicity carried by the 
charmed quark. Several remarks are in order: 

• Comparing the SU^'') (4) with SU^'') (4) predictions, 
one can see that there are 10-20% differences be- 
tween them only for the light-flavor contents, e.g. 
d-u (0.111 -> 0.141), d/u (0.71 -> 0.64), 2s/{u+d) 
(0.66 -> 0.75), and 2s/{u + d) (0.118 -> 0.133) 
[we note that both predictions are consistent with 
the existing data within errors] , and all other pre- 
dictions are almost no changes. It shows that the 
ry'" has minor or no effect on the spin and charm- 
related observables in the model. This is because 
first the 77'° does not have c, c components and sec- 
ond the q-f and (or and g^), q = u,d,s, ap- 
peared with equal probability in 77'°. It should be 
noted that the approximation C'^ — does not 
mean the probability of fluctuation q ^ q' + 77'" 
iC'^a) is smaller than that oi q ^ q' + D'^ (ccfl)- 
Actually, we have ec < C'^- Hence we use the ap- 
proximation C'^ = in SU'"''(4) version only for the 
practical reason explained here. 

• The probability a of the fluctuation u + 7r+ 
has been estimated in chiral field theory (see e.g. 

[§,13) 



where 



dze(A2-r(z))P„^rf+,+ (^), (49) 



i^u~*d+-K+\Z) — — Z 



32^2;.' 



gA — 0.75 is the dimensionless axial- vector cou- 
pling, /tt ~ 0.093 GeV the pion decay constant, 
{iTid) the constituent mass of the u (d) quark, 
A the ultraviolet cutoff, and m,^ is the pion mass. 
For A ~ 2.4 GeV, one obtains a ~ 0.142. For the 
strange quark fluctuation, u —^ s + K~^, similar 
calculation with some approximation leads to the 
probability ea ~ 0.062, which gives e ~ 0.44. We 
assume the same formula can be used for the charm 
fluctuation, u ^ c + D^, and obtain e^a ~ 0.004, 
which gives Cc — 0.03, which is consistent with 
ec ~ 0.06 ± 0.04 given in Eq.(45). Physically, 
the probability of splitting to the heavier mesons 
should be less than of splitting to the lighter ones. 
Hence the above estimation is quite reasonable. 

• The theoretical uncertainties shown in (46)- (48) 
and in Tables III, IV and V arise only from the 
uncertainty of in (45). If the observable does 
not depend on such as d — u, d/u, 2s/ {u + d), 
etc. (these quantities depend only on A, B, and 
G, and not on D, i.e. independent of e^), there is 
no uncertainty for them. This has been shown in 
Table III. Two special quantities Ac/c and As/s 
are also independent of ec as shown in Table IV. 
We put a star (*) mark on some 'data' in Tables 
III and IV to denote they arc model predictions or 
from theoretical analyses. 

• The SU(4) quark-meson model predicts an intrinsic 
charm component of the nucleon, 2c/ ^(g + g) ~ 
1%, which agrees with the predictions given in 

tand is also close to the those given in ||, |^ and 
p7[ . We note that the IC component is almost 
one order of magnitude smaller than the intrinsic 
strange component 2s/ + *?)• 

• Using the approach given in a previous work (see 
Eq. (3.6) in |2^), we can show that 

2 dxxc{x)/'S^ / dxx[q{x) + q{x)] 
Jo 

< 2c/^(g + g-), (51) 

where, as defined in this paper, c = dxc{x) and 

E(9 + 9), = Elod^i'iix) + q{x)]. The l.h.s. of 
Eq. (51) is the fraction of the total quark momen- 
tum carried by the charm and anticharm quarks. 
The prediction (46) implies that this fraction is less 
than 1%. Assuming the quark and antiquark share 
about one half of the nucleon momentum, then the 
charm and anticharm carry about 0.5% of the nu- 
cleon momentum or less. 



• From Table III, we have 



dt- 



A2 



[t-miy 



, (50) 



u + u 
- 0.53 



d + d 
0.37 : 



: s + s : c + c 
0.09 : 0.01. 



(52) 
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If we assume the quarks carry about 55% of the 
nucleon momentum, Eq. (52) imphes that the frac- 
tions of the nucleon momentum carried by u-, d-, 
S-, and c-flavors are approximately 29.2%, 20.3%, 
4.9%, and 0.6% respectively. They may compare 
with 31.4%, 17.8%, 4.3%, and 1.2% given by the 
DIS data at Q^=20 GeV^, where the gluons carry 
about 45% of the nucleon momentum. 

• The prediction of intrinsic charm polarization, 
Ac ~ -0.009 ± 0.006 is close to the resuh Ac = 
—0.020 ± 0.005 given in the instanton model 
This might be related to the strong suppression of 
non-spinflip contribution in both models. Our Ac 
in (47) is smaller in magnitude than that given in 
[5, yi (~ -0.3). However, the size of Ac given in 
[7| (~ —5 • 10~^) is even smaller. Hence further 
investigation in this quantity is needed. 

• The ratio Ac/ AY, as the function of ec is plotted in 
Fig.2. Taking ~ 0.06, one has Ac/AE ~ -0.021. 
This is consistent with the prediction given in 
but smaller than that given in . Combining with 
the fractions of the light quark helicities, we have 



Au/AE~ 2.171, Ad/AE ~ -0.988, 



As/AE ~ -0.162, Ac/AE ~ -0.021, 



(53) 



one can see that the u-quark helicity is positive 
(parallel to the nucleon spin) and about two times 
larger than the total quark helicity AE. The d- 
, S-, and c-helicities, however, are all negative 
[antiparallel to the nucleon spin), and their sizes 
are decrease as 



|Ad| : |As| : |Ac| ~ 1 : 10" 



10" 



(54) 



Compare to the intrinsic strange helicity As, the 
intrinsic charm helicity is one order of magnitude 
smaller. 

Since the hyperon /3-decay data are measured at low 
and and our model predictions are defined at the scale 
(0.2 GeV)^ < <(1.0 GeV)^, hence we may compare 
them with less ambiguity. However, many data listed in 
Tables HI and IV are coming from the DIS measurements 
at higher range. To make a meaningful comparison of 
model predictions with these data, we have to discuss the 



dependence of these observables. For spin observable, 
as we mentioned in section HI of that the model 
predictions, e.g. Au, Ad, As, etc. are compared with 
the (factorization) scheme-independent DIS observables 
a,(g2) = Ag- [as(g2)/27r]AG'(g2) {q = u,d,s) at the 
same scale [i.e. (0.2 GeV)^ < (^i q GeV)^], where 
AG{Q^) is the helicity of the gluon, and the axial charge 
o-qiQ^) is defined in the Adler-Bardeen scheme. Although 
aq{Q^) is independent of at the leading order and 
changes very slowly with at NLO, we still need to 
assume the perturbative QCD can be used down to the 
scale fj,^. The perturbative QCD evolution approach has 
been successfully used down to — 0-23 GeV^ Q, it 
is not clear, however, if the approach still hold below this 

Finally, in the quark meson fluctuation model, it is 
possible to include the contributions come from quark 
splittings to the vector mesons such as K* , p, etc. [1~ 
nonet in SU(3) and 1~ 15-plet in SU(4)]. However, it 
will change the formalism completely and is far beyond 
the goal of this paper. For simplicity and consistency, we 
only discuss the contributions of 0~ pseudoscalar mesons 
at this moment and defer the discussion of possible vector 
meson fluctuations to a later time. 

In summary, we have calculated the intrinsic charm 
contribution in the SU(4) quark meson model with sym- 
metry breaking. Despite the approximations and possi- 
ble theoretical uncertainties, the overall agreement be- 
tween the predictions and the existing data seems to be 
quite satisfactory considering the model is simple and 
has only a few parameters. The model also leads to 
many nevu predictions on observables explicitly related 
to the charm or anticharm. These observables are zero, 
e.g. 2c/ ^(g -I- g), Ac, etc. or indefinite, e.g. Ac/c and 
c^/c^, in the SU(3) description. We hope that these pre- 
dictions can be tested by the analyses of the DIS data on 
polarized and unpolarizcd charm productions in the near 
future. 
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Tables 



TABLE I: The probabilities Pq.|. (g^ | , | ) and Pg.^ (g^ | , g^ |). 



<? 










1 - (1 + e + ec + 2/)a/2 + (1 - Afa/2 


A^a/2 


&-a/2 




(1 + e + ec + 2/)a/2 + (1 - Afa/2 


a + A^a/2 


ea + B'^a/2 




A^a/2 


1 - (1 + e + + 2/)a/2 + (1 - Afa/2 


&-a/2 




a + A'^a/2 


(1 + e + + 2/)a/2 + (1 - Afa/2 


ea + B'^a/2 


St 


B^a/2 


B^a/2 


1 - (e + ^ + ec/2)a + (7^0/2 


Si 


ta + B^a/2 


ea + B^a/2 


(e + /. + ec/2)a + (7^0/2 


CT 


&al2 


D'^a/2 


£»2a/2 


Ci 


eca + b'^a/2 


e^a + D^a/2 


eca + I>^a/2 




(l~ Afa/2 


{l + A^)a/2 


[e + B^)a/2 


d],i 


(l + i2)„/2 


{I- Afa/2 


{e + B^)a/2 




{e + B'^)a/2 


{e + B^)a/2 






(ee + b'^)a/2 


(fe + 02)a/2 





TABLE II: Parameter sets for different symmetry-breaking schemes 



Model 


a 


e 


^12 

S S!7(3) 




SU<"'(4) 


0.143 


0.454 




0.06 ± 0.04 


SU('''(4) 


0.143 


0.454 


0.1 


0.06 


SU(3)|28| 


0.145 


0.460 


0.1 







TABLE III: Quark Flavor Observables. 



Quantity 


Data 


SU(")(4) 


SU('')(4) 


SU(3)pl 


NQM 


d — u 


0.110 ± 0.0181 
0.147 ± 0.0391 


35 
" 




0.111 


0.141 


0.143 





u/d 


[u{x)/d{x)]o.i<:c<o.2 = 0.67 ± O.OeH 
[u{x)/d{x)]^^o.i8 = 0.51 ± 0.06^ 


0.71 


0.64 


0.64 


— 


2s/{u + d) 


< 2xs{x) >/< x{u{x) + d{x)) > 
= 0.477 ±0.051|p8[ 


0.66 


0.75 


0.76 


- 


2c/(u + d) 




0.083 ± 0.055 


0.085 





- 


2s/{u + d) 


< 2xs{x) > / < x{u{x) + d{x)) > 
= 0.099 ± 0.009|p8[ 


0.118 


0.133 


0.136 





2c/(u + d) 




0.015 ±0.010 


0.015 












0.530 ± 0.004 


0.519 


0.523 


2/3 


{d + d)/Y.(q + q) 




0.368 ± 0.003 


0.370 


0.374 


1/3 


(s + -s)/j:{q + q) 


< 2xs{x) >/J2 < ^(<]i- 
= 0.076 ± 0.02^ 
= 0.10 ±0.06 
= 0.15 ±0.03 


- 

4C 


+ q{x)) > 

1 

1 

1 


0.090±0.001 


0.100 


0.103 





{c + c)/J2{q + q) 


0.03 |ll* 
0.02 ji* 
0.01 p]* 
0.009 J 1* 
0.005 [§|* 
<0.004 ||l7[* 


0.011 ±0.008 


0.011 








El/El 


J2 < xq{x) >/J2 < xq{x) > 
= 0.245 ± 0.005@ 


0.230 ± 0.004 


0.233 


0.231 






TABLE IV: Quark Spin Observables 



Quantity 


Data 


SU("'(4) 


SU(''H4) 


SU(3)|25] 


NQM 


Au 
Ad 
As 


0.85±0.04| 
-0.41±0.04 
-0.07±0.04 


4l[ 

B 


0.871 ± 0.009 
-0.397 ± 0.002 
-0.065 ± 0.000 


0.859 
-0.393 
-0.065 


0.863 
-0.397 
-0.067 


4/3 
-1/3 



Ac 


-0.020 ±0.004 @* 
-0.3 1* 

-5 -10-* ir 


-0.009 ±0.006 


-0.009 








AE/2 


0.19±0.06pll 


0.200 ± 0.006 


0.196 


0.200 


1/2 


Au, Ad 


-0.02 ±0.11(42 














As, Ac 
















Au/AE 
Ad/AE 
As/AE 




2.171 ±0.043 
-0.988 ± 0.024 
-0.162 ±0.005 


2.192 
-1.004 
-0.166 


2.162 
-0.994 
-0.167 


4/3 
-1/3 



Ac/AE 


-0.08 ±0.01 §* 
-0.033 1^* 


0.021 ± 0.014 


-0.022 








Au/u 




0.383 ± 0.003 


0.381 


0.383 


2/3 


Ad/d 




-0.287 ± 0.001 


-0.283 


-0.284 


-1/3 


As/s 




-3/10 


-0.269 


-0.269 




Ac/c 




-16/51 


-16/51 










2.241 ± 0.012 


2.231 


2.241 


5 


d^/dl 




0.554 ±0.001 


0.559 


0.558 


1/2 






7/13 


0.576 


0.576 




ct/c; 




35/67 


35/67 
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0.136 zfco.oiepi 


0.143 ±0.002 


0.141 


0.142 


5/18 




-0.041 ± 0.007§|] 


-0.042 ± 0.001 


-0.043 


-0.042 





as 


1.2670±0.0035P] 


1.268 ±0.010 


1.252 


1.260 


5/3 


ag 


0.579± 0.0253 


0.605 ± 0.006 


0.595 


0.600 


1 



TABLE V: F, D and Ga/Gv ratios of the hyperon beta decays, [a]: This is the best fit to four measured (Ga/Gv) (x^ — 1-96) 



under tlie constraint F±D=1.267. 



Quantity 


Data 


SU("'(4) 


SU<'''(4) 


SU(3)H 


NQM 


F±D 


1.267[a] 


1.268 ±0.010 


1.252 


1.260 


5/3 


F 


0.463[a] 


0.468 ± 0.004 


0.462 


0.465 


2/3 


D 


0.804[a] 


0.800 ± 0.006 


0.790 


0.796 


1 


F/D 


0.576[a] 


0.585 ± 0.007 


0.585 


0.585 


2/3 


{GA/Gv)n^p 


1.2670 ± 0.00350 


1.268 ±0.010 


1.252 


1.260 


5/3 


{Ga/Gv)k^p 


0.718 ±O.O150r 


0.735 ± 0.007 


0.725 


0.730 


1 


{GA/Gv)s-~^n 


-0.340 ± 0.0171 




-0.332 ± 0.002 


-0.328 


-0.330 


-1/3 


{Ga/Gv)s- -,A 


0.25±0.05[p3 


1 


0.202 ± 0.002 


0.198 


0.200 


1/3 



